arXiv:1509.05628v2 [astro-ph.HE] 30 Dec 2015 


MNRAS 000, 1-9 (2015) 


Preprint 31 December 2015 


Compiled using MNRAS style file v3.0 


Positive correlation between the cyclotron line energy and 
luminosity in sub-critical X-ray pulsars: Doppler effect in 
the accretion channel 

Alexander A. Mushtukov,^’^* Sergey S. Tsygankov,^ Alexander V. Serber,^ 
Valery F. Suleimanov and Juri Poutanen^’® 

^ Tuorla observatory, Department of Physics and Astronomy, University of Turku, Vdisaldntie 20, FI-21500 Piikkid, Finland 
^Pulkovo Observatory, Russian Academy of Sciences, Saint Petersburg 196140, Russia 
^Nizhny Novgorod Planetarium, Revolutsionnaya st., 20, Nizhny Novgorod 603002, Russia 
Institut fiir Astronomic und Astrophysik, Universitdt Tubingen, Sand 1, D-72076 Tubingen, Germany 
^ Kazan (Volga region) Federal University, Kremlevskaja str., 18, Kazan 420008, Russia 

^Nordita, KTH Royal Institute of Technology and Stockholm University, Roslagstullsbacken 23, SE-10691 Stockholm, Sweden 


Accepted 2015 September 18. Received 2015 September 17; in original form 2015 July 17 


ABSTRACT 

Cyclotron resonance scattering features observed in the spectra of some X-ray pul¬ 
sars show signihcant changes of the line centroid energy with the pulsar luminosity. 
Whereas for bright sources above the so called critical luminosity these variations are 
established to be connected with the appearance of the high accretion column above 
the neutron star surface, at low, sub-critical luminosities the nature of the variations 
(but with the opposite sign) has not been discussed widely. We argue here that the 
cyclotron line is formed when the radiation from a hotspot propagates through the 
plasma falling with a mildly relativistic velocity onto the neutron star surface. The 
position of the cyclotron resonance is determined by the Doppler effect. The change of 
the cyclotron line position in the spectrum with luminosity is caused by variations of 
the velocity profile in the line-forming region affected by the radiation pressure force. 
The presented model has several characteristic features: (i) the line centroid energy 
is positively correlated with the luminosity; (ii) the line width is positively correlated 
with the luminosity as well; (iii) the position and the width of the cyclotron absorp¬ 
tion line are variable over the pulse phase; (iv) the line has a more complicated shape 
than widely used Lorentzian or Gaussian profiles; (v) the phase-resolved cyclotron line 
centroid energy and the width are negatively and positively correlated with the pulse 
intensity, respectively. The predictions of the proposed theory are compared with the 
variations of the cyclotron line parameters in the X-ray pulsar GX 304-1 over a wide 
range of sub-critical luminosities as seen by the INTEGRAL observatory. 
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1 INTRODUCTION 

X-ray pulsars (XRPs) are neutron stars (NSs) in binary sys¬ 
tems accreting matter usually from a massive optical com¬ 
panion. Because of the strong magnetic field B > 10^^ G) 
the accreting plasma is channelled toward the NS magnetic 
poles, where the accretion flow heats the surface and most of 
the gravitational energy is released via hard X-ray emission. 
The magnetic field modihes the observed X-ray spectrum 
often manifesting itself as the line-like absorption features, 
the so-called cyclotron lines. Such cyclotron resonance scat- 

* E-mail: al.mushtukov@gmail.com (AAM) 


tering features, sometimes also with harmonics, are observed 
in the spectra of about two dozens of XRPs (Coburn et al. 
2002; Filippova et al. 2005; Walter et al. 2015). 

In some cases, the luminosity-related changes of the 
line energy are observed, suggesting that the configura¬ 
tion of the line-forming region depends on the accretion 
rate. The line energy has been reported to be positively 
(for relatively low-luminous sources, L < lO^^ergs”^; 
see Staubert et al. 2007; Klochkov et al. 2012) and neg¬ 
atively correlated with luminosity (for relatively high- 
luminous sources, L > 10®'^ erg s“^; see Mihara et al. 1998; 
Tsygankov et al. 2006, 2007, 2010), as well as uncorrelated 
with it (Caballero et al. 2013). 
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Figure 1. The schematic presentation of the dependence of the cyclotron line energy on the velocity profile in the line-forming region. 
The radiation pressure affects the velocity profile when the luminosity is approaching its critical value: the higher the luminosity, the 
lower the velocity in the vicinity of NS surface. The lower the electron velocity, the lower the redshift and the higher the resonant energy 
at a given height. As a result, a positive correlation between the line centroid energy and the luminosity is expected. Panel (a) corresponds 
to the low-luminosity case when the gas is in a free-fall. The observed cyclotron line is redshifted relative to the rest-frame value. Panel 
(b) is for the higher luminosity when the gas is decelerated close to the surface and the line shifts closer to the rest-frame position. 


The difference between positively and negatively cor¬ 
related cases is explained by the absence or presence of 
the accretion column above the stellar surface, which be¬ 
gins to grow as soon as the accretion luminosity reaches 
its critical value L* ~ 10®'^ erg s“^ (Basko & Sunyaev 1976; 
Mushtukov et al. 2015b, c). The negative correlation in a 
high-luminosity case was discussed in a number of works: 
Poutanen et al. (2013) proposed that the line is produced in 
the spectrum reflected from the NS surface and the anti¬ 
correlation with luminosity is explained by variations of 
the NS area illuminated by a growing accretion column. 
Nishimura (2014) suggested that the growth of the column 
is associated with the decreasing magnetic field that leads 
to a shift in the cyclotron line position. This requires, how¬ 
ever, that the column is only a few hundred meters high, 
contradicting the theoretical estimates giving a much taller 
column. The positive correlation in a low-luminous case has 
been explained by variations of the atmosphere height above 
the NS surface because of the changing ram pressure of 
the infalling material (Staubert et al. 2007). This model is 
based, however, on an ad hoc assumption that the stop¬ 
ping depth of accreting protons (identified there with the 
atmosphere scale-height) is -^lOOm instead of a more real¬ 
istic value of 1 m. Alternatively, Nishimura (2014) proposed 
that changes in the emission pattern towards more fan-like 
diagram at higher luminosity affect the position of the cy¬ 
clotron line via the Doppler effect. There it was assumed 
that the velocity profile is not affected. This is, however, 
a gross simplification, because the accretion velocity in the 
vicinity of the NS surface is expected to change from the 
free-fall value of ~ c/2 to 0 with the grow of the luminosity 
to the critical value L*. Additionally, this model seems to 
contradict the positive correlation of the line width with lu¬ 
minosity (Klochkov et al. 2012), as the thermal width of the 
line is expected to be much smaller for directions across the 
magnetic field because of a weak transverse Doppler effect 
(Mushtukov et al. 2015a). 


In this work we will focus only on the low-luminosity 
case (sub-critical XRPs). In this case, bulk of the radia¬ 
tion comes from the hotspot heated by the infalling plasma 
(Zel’dovich & Shakura 1969; Nelson et al. 1993). This radia¬ 
tion, however, has to pass through the material falling inside 
the accretion channel. The optical thickness of the channel 
at the resonant energies is much higher than unity in the 
considered luminosity range, L ~ 10^® fi- 10®^ erg s“^. The 
photons are resonantly scattered there, leading to the forma¬ 
tion of a cyclotron absorption-like feature in the spectrum. 
Plasma in the accretion flow moves with a significant veloc¬ 
ity towards the NS surface, and the scattering feature must 
be redshifted due to the Doppler effect. An observer thus 
sees the line at a lower energy than the cyclotron one asso¬ 
ciated with the magnetic field at the pole. The amplitude 
of the shift depends on the electron velocity and the an¬ 
gle between the photon momentum and the accretion flow 
velocity. As a result, the observed position and the shape 
of the cyclotron scattering feature depend on the velocity 
profile in the accretion channel (see Fig. 1). Particularly, 
it is expected that the lower the electron velocity in the 
region, the lower the redshift and the higher the line cen¬ 
troid energy. In this case, it is also closer to the actual 
cyclotron energy E^yc — 11.6i3i2 keV, which is defined by 
the surface magnetic field strength B 12 = The 

radiation pressure affects the velocity profile when the lu¬ 
minosity is approaching its critical value (Basko & Sunyaev 
1976; Mushtukov et al. 2015b): the higher the luminosity, 
the lower the velocity in the vicinity of NS surface. There¬ 
fore, we expect a positive correlation between the line cen¬ 
troid energy and the XRP luminosity. 

In this paper, we construct a simple model, which de¬ 
scribes variations of the opacity near the cyclotron energy, 
and compare the predictions of our model with the observa¬ 
tions obtained by the INTEGRAL space observatory of the 
XRP GX 304-1 during its Type II outburst in the beginning 
of 2011. 
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Figure 2. The dependence of (a) the dimensionless velocity and 
(b) the corresponding resonant scattering energy in the laboratory 
reference frame on the distance from the NS surface measured 
in units of the thickness of the accretion channel. The velocity 
profile is given by equations (6) and (7) and the cyclotron energy 
is computed from equation (1). Different curves correspond to 
various final electron velocities: /3o = 0 (solid), /3o = 0.2 (dashed), 
/3o = 0.4 (dotted). Here, the free-fall velocity is /Sg = 0.5 and the 
observing angle is 6 = 0. 


2 THE MODEL 

The basic process which determines the interaction between 
the radiation and matter in the vicinity of a NS surface is 
Compton scattering. In case of strong magnetic field the 
scattering has a number of special features. Its cross section 
strongly depends on the B-field strength, photon energy, 
polarization and momentum. Electron transitions between 
Landau levels cause the resonance scattering at some pho¬ 
ton energies, where the cross section can exceed the Thom¬ 
son scattering cross section by several orders of magnitude. 
The resonant scattering leads to appearance of cyclotron ab¬ 
sorption features in the spectra of XRPs. The exact shape 
and position of the features in the spectra are determined 
by the structure of the line-forming region. 

The exact resonant energy depends on the velocity of 
the electrons because of the effects of special relativity. In 
a particular case of cold electrons in bulk motion with the 
velocity /3 = u/c towards the NS surface, the position of the 
fundamental is (Daugherty & Harding 1986; Meszaros 1992; 


Harding & Lai 2006) 

_ CbC-P cos 9')]^ + 2b si^ - 7(1 - /I cos 6 )') 
meC^ sin^ 6 ' 

( 1 ) 

where 7 = (1 — /3^) is the Lorentz factor, 9 ' is the angle 
between the electron and photon momentum, b = BjBci 
is the local magnetic field strength in units of the critical 
magnetic field Bcr = /(eh) = 4.413 x 10^® G. For 6 -C 1 

or 0 ' <C 1 , the expression for the fundamental in the observer 
(laboratory) reference frames is reduced to 


_ b 

rrieC^ 7 ( 1 -I-/3cos0) ’ 


( 2 ) 


where 6 — n —9' is the angle between the line of sight and the 
magnetic axis (normal to the NS surface). We see that the 
resonant energy in the laboratory reference frame is always 
redshifted and the redshift depends on the plasma velocity 
and the observing angle. 

The velocity profile of the infalling material in the vicin¬ 
ity of the NS surface is determined by the radiation flux 
F and the effective cross-section a, because the accreting 
matter is interacting with radiation and is decelerating to 
some final velocity jSo = vq/c at the NS surface. Then 
it looses its momentum entirely due to Coulomb collisions 
(Zel’dovich & Shakura 1969). The velocity just above the 
NS surface /3o and the velocity profile in the accretion chan¬ 
nel depend on the radiation pressure force Fa/c and there¬ 
fore on the luminosity of an XRP. Obviously, the velocity fio 
tends to be smaller for a higher luminosity and it turns to 
zero when the XRP luminosity reaches its critical value 


L* = 4 X 10 ®« ^ 


I 


Rs \2 X 10®cm 


(Tt -1 

— erg s , 
a 


(3) 


where m = M/Mq is the NS mass in units of the solar 
mass, R = 10®Re cm is the NS radius, I is the length of 
the accretion channel footprint, which has a shape of thin 
annual arc, and ctt is the Thomson cross-section. The crit¬ 
ical luminosity L* strongly depends on the R-field strength 
(Basko & Sunyaev 1976; Mushtukov et al. 2015b). At higher 
luminosity the accretion column begins to grow above the 
surface. 

If the luminosity L ^ L*, the gravitational force is much 
lower than the radiation pressure force, which defines the 
deceleration of matter. In order to obtain the velocity profile, 
we need to account for the dependence of the radiation flux 
on the height h above the NS surface. In case of a filled 
accretion funnel, the dependence is 


where Fq is the flux at the surface and d here is the ra¬ 
dius of the hot spot. Then the velocity changes with height 
according to the equation 

^ _ 2 GF{h) 

Ah m-pC® ’ 

where rrip is the proton mass. Thus we get (see Fig. 2) 

1/2 


(5) 


/3(h) = 
where 
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If /Jo = 0, then _Fb becomes the critical flux F* corresponding 
to the critical luminosity L* = 2tv<AF*. It is also obvious 
that the luminosity corresponding to a given velocity /3o at 
the bottom of the accretion channel is 

L(/3o) L‘ (^1 - . (8) 

If the accreting matter is confined to a narrow wall 
of the magnetic funnel, the flux dependence on the height 
is more complicated. The flux drops as (1 + h? 
when h -C l/{‘2'!r) and as (1 + when h > 

l/{2n), resulting in a slightly different velocity profile. 
However, the qualitative behaviour is the same (see also 
Lyubarskii & Syunyaev 1982; Becker 1998) and further in 
this paper we use equation (6) for calculation of the velocity 
profile. 

Radiation from the hot NS surface has to pass through 
the layer of moving plasma and is resonantly scattered there 
by electrons. The resonant energy depends on the local elec¬ 
tron velocity (see Fig. 2). As a result, it is expected that the 
position of the absorption line feature in the spectrum and 
its shape are determined by the plasma velocity profile near 
the NS surface. Higher luminosity would lead to lower veloc¬ 
ity shifting the resonance to higher energy in the observer 
reference frame. 


2.1 The opacity of moving plasma near the 
resonance energy 

In order to dehne main characteristics of the cyclotron line 
and its variability with changes of the accretion luminosity, 
we have to calculate the opacity of the line-forming region 
as a function of photon energy E. The opacity in a given 
direction (defined by the zenith angle 0 and the azimuthal 
angle <j)) for photons of energy E emitted at the distance x 
from the accretion channel edge is 

T{E,e,<l>,x)^ I ^a{E,e,y)nAy), (9) 
0 

where y*{9,(j)) ~ x /(tan 6 cos 4>) for the case of J d (in our 
case l/d ~ 160L37 = L/lO^^erg s“^, 
see Mushtukov et al. 2015c), cr{E, 9, y) is the scattering cross 
section at height y above the stellar surface, Ue(i/) is the 
local electron number density. The scattering cross section 
cr{E,9,y) depends on the local temperature T{y) as well. 
We use the scattering cross-section computed using sec¬ 
ond order perturbation theory of quantum-electrodynamics 
(Daugherty & Harding 1986; Harding & Daugherty 1991; 
Mushtukov et al. 2012, 2015a). 

The resonant scattering cross-section of the X-mode 
photons is higher than the resonant scattering cross-section 
of the 0-mode photons (Canuto et al. 1971; Nagel 1981; 
Mushtukov et al. 2015a). As a result, the emission around 
the cyclotron line should be dominated by the O-mode pho¬ 
tons. Therefore, we use the O-mode cross-sections in the fol¬ 
lowing calculations. We also assume that the electron tem¬ 
perature of the accreting gas is constant and we fixed it 
at 3keV, a value close to the Compton temperature of ra¬ 
diation. The results are not affected by the specific choice 
of that temperature. Then the average value of the opacity 


over the coordinate x inside the accretion channel and the 
azimuthal angle is 

d 7r/2 

t{E,9)~-^ J dx J d4>T{E, 9, (j>,x)exp{-TT{9,(j>,x)), 

0 0 

( 10 ) 

where rT{9,(j>,x) is the Thomson optical thickness at the 
direction given by angles 9 and (j), which affects the resid¬ 
ual continuum intensity at a given direction from given 
point: I{9,(j>) « l{9,(j), x) exp{—TT{9,(fi,x)). If the accre¬ 
tion disc is interrupted in the C-zone (Shakura & Sunyaev 
1973; Suleimanov et al. 2007), where the gas pressure and 
the Kramer opacity dominate, the accretion channel thick¬ 
ness can be estimated as follows: 

d « 6 X 10^ cm, (11) 


with A = 0.5 being a commonly used value (Ghosh & Lamb 
1978). 

In order to compare with observations, we then compute 
the opacity averaged over the pulsar phase y>: 


^ £ t{E, 9)\ cos9\dip 

J,Acos9\dp 


( 12 ) 


where cos 9 = cos i cos a + sin i sin a cos (p, i is the angle be¬ 
tween the line of sight and the NS rotational axis, a is the 
magnetic inclination and (p is the phase angle. Here we as¬ 
sumed a blackbody-like emission pattern with the emitted 
flux proportional to the cosine of the hotspot inclination to 
the line of sight cos 9. 

Assuming dipole magnetic field structure, the field 
strength in the vicinity of stellar surface depends on height 

h as follows: B{h) = Bq , where Bo is the surface 

magnetic field strength. From the mass conservation law 
M/(2Sd) — mpne{h)v{h) we can get the local electron num¬ 
ber density: 


ne(h) « 2x10^° ^^^^'^^^^ « 2x10'® L%^B^^^'^p-\h) cm-^ 

(13) 

where P{h) = v{h)/c is dimensionless velocity and is the 
hotspot area: 

« 3 X 10® A-^/« cm®. (14) 


Thus for a given luminosity, we compute the velocity pro¬ 
file in the line-forming region, which affects the opacity, and 
then for given parameters i and a we can compute the av¬ 
erage opacity using equation (12). 


2.2 Opacity profiles 

For our calculations we have taken a NS of mass M = 
1.4 Mq, radius R = 10 km and surface dipole magnetic field 
strength at the poles Bq = 5x lO'® G. The opacity averaged 
over the orientation of line-forming region (see eq. (12)) as a 
function of photon energy is presented in Fig. 3, where dif¬ 
ferent lines correspond to different final velocity on the NS 
surface. The peak of the opacity, whose complicated shape 
is defined by the velocity distribution in the line-forming re¬ 
gion, shifts to higher energies when the surface velocity of 
the accretion flow decreases. The lower the surface velocity, 
the wider the opacity peak produced by resonant scattering. 
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Figure 3. The averaged opacity profile for various final velocities 
/3o at the NS surface. The profiles are computed for the fixed 
accretion luminosity L = 10^^ erg s“^. We see that the lower the 
final velocity, the higher the centroid energy and the wider the 
feature. Here Bq = 5 X 10 ^^ G, i = 60 °, a = 20 °. 

The opacity profile determines the shape of the scattering 
feature in the XRP spectrum. Therefore, taking into account 
the connection between the accretion luminosity and the ve¬ 
locity j3o we find that the cyclotron absorption feature cen¬ 
troid energy and its width should be positively correlated 
with the luminosity for sub-critical XRPs: the higher the lu¬ 
minosity, the higher the line centroid energy, the wider the 
cyclotron line in the spectrum. The increase of the centroid 
energy is caused by the appearance of a layer with a lower 
velocity, where the redshift of the cyclotron energy is lower. 
This fact explains also the increase of the line width: it grows 
on its blue side because of the layers of lower velocity. We 
also note that the shape of the line is asymmetric and can 
hardly be approximated by the widely used Lorentzian or 
Gaussian profiles. 

The behaviour of the line can be well illustrated by the 
ratio of opacities calculated for different surface velocity val¬ 
ues (see Fig. 4). If one takes the ratio of the low-luminosity 
opacity (high /3o) to the high-luminosity opacity (low /3o), 
then the shape of the ratio is quite conservative and does 
not depend on the angles, which the opacity is averaged over: 
the ratio reaches a local maximum and then has a deep local 
minimum (see Fig. 4b). The same behaviour is expected for 
the ratio of higher luminosity XRP spectrum to the lower 
luminosity spectrum. Because we do not have a model for 
the spectrum formation, the analysis of spectral ratios can 
be suggested as a good test for the model. 

It is also expected that the cyclotron line absorption 
feature has to show variability during the pulse period. The 
Doppler shift of the resonant energy is lower for larger an¬ 
gles between the accretion flow velocity and the line of 
sight (see eq. (2)). Therefore, it is expected that the phase- 
resolved line centroid energy would be anti-correlated with 
the pulse intensity for sub-critical XRP with a pencil beam 
emission diagram (Gnedin & Sunyaev 1973). Because the 
resonance peak in the scattering cross-section tends to be 
wider for the photons propagating along the B-held lines 
(Harding & Daugherty 1991; Mushtukov et al. 2015a), it is 
expected that the phase-resolved cyclotron line width would 



E [keV] 

Figure 4. (a) The opacity profile averaged over the pulse period 
for various magnetic inclination a (see eq. (12)). Black lines are 
for i = 60° and 0 = 5°, while red lines are for i = 60° and a = 
20°. Solid and dotted lines are given for do = 0-1 and do = 0-4, 
respectively, (b) The ratio of the opacity profiles ti = r(do = 0.4), 
T 2 = 'r(do = 0.1) corresponding to the different velocities at the 
NS surface. Here B = 5 X 10^^ G. 

be positively correlated with the pulse intensity. The radia¬ 
tion beam pattern also affects the phase-averaged spectrum 
and its changes with the luminosity can affect the cyclotron 
line energy (Nishimura 2014). In our calculation the radi¬ 
ation beam pattern is taken into account in equation (10), 
where the opacity is weighted with the continuum intensity 
escaping in different direction. 

3 COMPARISON WITH OBSERVATIONS 
3.1 Observational data 

The most reliable positive correlation between the cy¬ 
clotron line energy and the luminosity has been observed 
in the spectra of XRP GX 304-1 (Yamamoto et al. 2011; 
Klochkov et al. 2012). Therefore, we use this source as a 
case study to verify the predictions of our model. 

GX 304-1 belonging to the subclass of XRPs with Re¬ 
type optical companions (Mason et al. 1978) has a relatively 
long pulse period of 272 s (McGlintock et al. 1977). It is lo¬ 
cated at a distance of about 2.4 kpc (Parkes et al. 1980). 
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The cyclotron absorption line at 54 keV in the source energy 
spectrum was discovered only recently (Yamamoto et al. 
2011) using Suzaku and RXTE data obtained during Type I 
outbursts recurrently observed from this system every 
132.5 d. 

In our current work, we reanalyse the data obtained by 
the INTEGRAL observatory (Winkler et al. 2003) during 
another Type II outburst occurred in the beginning of 2011 
(Yamamoto et al. 2012). The data sample is identical to the 
one utilized by Klochkov et al. (2012). The total duration of 
the outburst was about one month. It was evenly covered 
by eight INTEGRAL observations (during spacecraft revo¬ 
lutions 1131-1138). Both rising and declining parts as well 
as the maximum of the outburst were monitored. 

The data from the IBIS/ISGRI detector (Lebrun et al. 
2003) were screened and reduced following the procedures 
described in Krivonos et al. (2010) and Churazov et al. 
(2014). Energy calibration was provided by the Offline Sci¬ 
entific Analysis (OSA) version 10.1 (Caballero et al. 2012). 
Additional gain correction based on the position of the tung¬ 
sten fluorescent line (at 58.2 keV) was applied, keeping the 
uncertainty of energy reconstruction to less than 0.1 keV. 
Data from the JEM-X monitor (Lund et al. 2003) were pro¬ 
cessed with the standard procedures from OSA 10.1. 

Observed fluxes from GX 304-1 range between about 
0.1 and 1 Crab, permitting us to trace the evolution of 
the energy spectrum in a wide energy range (5-100 keV) 
over very different luminosities of the source. Three spec¬ 
tra obtained at different states are shown in Fig. 5(a). 
Particularly, spectra obtained during the lowest-luminosity 
state at the rising phase, peak phase, and lowest-luminosity 
state at the declining phase of the outburst are denoted by 
Rev. 1131, Rev. 1134, and Rev. 1138, respectively. JEM-X 
data are used in the range 5-30 keV and shown by green 
and red circles (for JEM-Xl and JEM-X2 modules, corre¬ 
spondingly), IBIS/ISGRI data cover the energy range from 
25 to 100 keV (shown by filled black circles). The solid 
curves represent the best-fitting models that consist of a 
power law with an exponential cutoff modified by the Gaus¬ 
sian line at 6.4 keV and the cyclotron absorption feature 
in the form of Lorentzian optical depth profile (cyclabs 
model in XSPEC; Mihara et al. 1990). The final form of 
the fitting model in XSPEC package can be expressed as 
(gau+powerlaw*highecut)*cyclabs. 

Three lower panels of Fig. 5 show the ratio of spectra ob¬ 
tained at different orbits (see Table 1). We see that the spec¬ 
tra during brighter states are harder. Such a behaviour can 
be accounted for by the increasing Compton j/-parameter 
with increasing mass accretion rate (Postnov et al. 2015). 
The most important feature in all spectral ratios is a wave¬ 
like structure at energies between 40 and 80 keV. We see 
an excess and depression relative to the general trend below 
and above 55 keV, respectively. 

Our results confirm conclusions made by 
Klochkov et al. (2012) about the positive correlation 
of the cyclotron energy with flux in GX 304-1. Dependence 
of the cyclotron line centre on the observed flux is shown 
in Fig. 6(a) for both shapes of the optical depth profile, 
Lorentzian (cyclabs model; red points) and Gaussian 
(gabs model; blue points). The corresponding dependences 
of the line width on flux are shown in Fig. 6(b). We see 
that both the line energy and its width correlate with the 



6 8 10 20 40 60 80 100 

Energy, keV 


Figure 5. (a) Energy spectra of GX 304-1 obtained by main in¬ 
struments of the INTEGRAL observatory during different space¬ 
craft revolutions. Data in the energy range 5—30 keV are from 
JEM-X monitor (red and green points for JEM-Xl and JEM-X2 
modules), between 25 and 100 keV — from the IBIS telescope. 
(b)-(d) Ratio of spectra obtained during different revolutions. 


flux. However, it is very important to keep in mind that the 
opacity profile differs significantly from a simple Gaussian 
or Lorentzian (see Sect. 2.2) at any intensity levels. These 
deviations leads to a difficulty with interpretation of the 
cyclotron energy obtained directly from the fits using such 
simple phenomenological models. 


3.2 Observations versus theory 


For given geometrical parameters i, a (see Section 2.1) and 
accretion luminosity L (or corresponding final velocity Pq) 
we get the averaged opacity as a function of energy (see 
Fig. 3). Then the cyclotron line centroid energy can be esti¬ 
mated as 


Ecyc (L) 


f ET(E,L)dE 
jT{E,L)dE ■ 


(15) 


For the particular set of parameters {i = 70°, a = 10°, L* = 
2.7 X 10®^ergs“^ and surface cyclotron energy Ecyc,o = 
66.5 keV) the theoretical dependence of the line centroid en¬ 
ergy on luminosity is shown in Fig. 6(a). It is qualitatively 
consistent with the observed variations of the cyclotron line 
centroid with luminosity. 

Because there is no model of the line formation, we can¬ 
not predict well the shape of the absorption feature and its 
width. The line width in the observed spectrum can be dif¬ 
ferent from the width of the resonance in optical thickness. 
We have measured the full width of the theoretical t{E) de¬ 
pendence at the level of 10 and 25 per cent of the maximal r 
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Table 1. Cyclotron line parameters for the INTEGRAL observations analysed in this work. 


Revolution 

T ® 


CYCLABS 



GABS 


number 

(10^^ erg s-l) 

Ecyc 

f^cyc 

Tcyc 

Ecyc 

f^cyc 

Tcyc 



(keV) 

(keV) 


(keV) 

(keV) 


1131 

0.50 


9-071.5 

U.DO_q 07 

56.sill 

7 k+0-9 
'•^-0.8 

10.411° 

1132 

1.12 

55.9t°j 

11.811° 

0-74l°;°® 

59.27°;° 

A 1+0.7 
^•-^-0.7 

14.3111 

1134 

2.20 

56.11°;® 

14.8ll» 

0-68l°;°° 

61.57°;? 

11.711? 

18.81?;? 

1136 

0.96 

52.6l°;° 

7-7713 

0-53l°;°® 

54.51°;? 

6.811? 

7.5l?-? 

1137 

0.47 

53.ilg;? 

5 1+1-5 

0.48l°:°« 

54.3l°;° 

5 . 911 ? 

5.5l?l 

1138 

0.22 

49.6li;° 

o a+2.5 

0.39l°-f° 

50.611° 

4 . 811 ? 

3.8l?;? 


“ Luminosity in the 3—100 keV energy band, assuming source distance of 2.4 kpc. 



I n a37 -It 

L [10 erg s ] 




T riA^^ -In 

L [10 erg s ] 


Figure 6. (a) The cyclotron line energy and (b) the cyclotron 
line width dependencies on the 3-100 keV luminosity of GX 304- 
1. Red squares show results obtained using the Lorentzian profile 
of the optical depth (cyclabs model); blue circles - using the 
Gaussian profile (GABS model) for the cyclotron absorption fea¬ 
ture. Luminosities are calculated assuming the distance to the 
source of 2.4 kpc. Theoretical relations are given by black solid 
lines for the cyclotron energy at the NS surface for parameters 
= 66.5keV, i = 70°, a = 10°, L* = 2.7 x lO^’’erg s”!. 
A pure 0-mode was considered. The width of the resonance is 
measured at the level of 10 and 25 per cent of the maximum 
t{E). 


Figure 7. Ratio of spectra of GX 304-1 obtained in revolutions 
1134 to 1138 corrected for different hardness ratio is given by blue 
circles (see also Fig. 5) and the ratio of the optical thicknesses for 
the case of pure O-mode polarization (solid line) for parameters 
E^ycQ = 66.5keV, C = 0.55, i = 70°, a = 10°. 


as a function of the accretion luminosity (see black solid lines 
in Fig. 6b). We see that the model predicts a positive corre¬ 
lation, which is consistent qualitatively with the data. More¬ 
over, we see that the observed ratios of the spectra obtained 
in different intensity states show very peculiar behaviour at 
the energies close to the cyclotron energy: the ratios reach 
local maximum below the cyclotron energy, then decrease 
and reach local minimum above the cyclotron energy (see 
Fig. 5). The same behaviour is seen in the ratios of the opac¬ 
ity profiles (see Fig. 4(b) and Fig. 7) calculated for different 
velocity gradients and, therefore, different accretion lumi¬ 
nosity. This observation is consistent with our model that 
the cyclotron line forms by scattering of the hard radiation 
from the hotspots off the electrons in the accretion flow and 
that the behaviour of the cyclotron line is determined by the 
changes of the velocity profile in the line-forming region. 

In order to estimate the depth of the absorption fea¬ 
ture, a detailed model of the line formation region is re¬ 
quired. Here instead we make base our analysis on simple at¬ 
mosphere models. According to the Schwarzschild-Schuster 
model (Mihalas 1978; Sobolev 1985), where the continuum 
and the absorption features form separately, the flux in the 


MNRAS 000, 1-9 (2015) 





























8 A. A. Mushtukov et al. 


vicinity of the scattering feature can be estimated as follows 

Fe cc (16) 

where te is the optical thickness of the line-forming re¬ 
gion for a given photon energy E. The Thomson optical 
thickness of the line-forming region is at L ~ L*, and 
therefore the resonant optical thickness of the region is 
much larger than nnity. As a result, the absorption fea¬ 
tures obtained with the Schwarzschild-Schuster model are 
extremely deep, which is not the case according to the ob¬ 
servational data (see Fig. 5(a) and Table 1). This discrep¬ 
ancy might be caused by the increasing flux towards the 
photosphere which is related to the energy-release inside 
the line-forming region where the accretion flow deceler¬ 
ates (Lyubarskii & Syunyaev 1982). Unfortunately, the ex¬ 
act energy deposition rate can be obtained only from a self- 
consistent model that includes solution of coupled hydrody¬ 
namic and radiation transfer equations. The photon redistri¬ 
bution over the energy can be also important. According to 
the Eddington model, where the regions of continuum and 
absorption feature formation are not separated, the absorp¬ 
tion feature depth in a case of monochromatic scattering is 
Te = FE/Fcont,E ~ ('Tcont.B/'TB)^/^ (Mihalas 1978; Sobolev 
1985), where Fe and Fcont.B are fluxes at a given energy 
E with and without the absorption feature, and Tcont.B is 
the optical thickness without resonant scattering. It gives 
r_B ~ 0.03 in our case. But if one takes into account photon 
redistribution over the energy, then te ^ {aT^ont,E 
(Sobolev 1985), where a = Trad/{2ujBVth/c) is the natu¬ 
ral line width parameter. Trad is the cyclotron decay rate, 
ojb is the cyclotron frequency and Uth = (2kBT. 
Frad/cus ~ 10~® for B = 6x 10^^ G (Herold et al. 1982) and 
therefore a ~ 10~^. As a result, the expected absorption 
feature depth te ~ 0.1. This is not far from the observed 
value (see Fig. 5). 

In order to compare theoretical and observational wave¬ 
like structures, we use the ratio of spectra obtained at revolu¬ 
tion numbers 1134 and 1138, when the accretion luminosities 
differ by a factor of 10, and correct the result for the differ¬ 
ence in the spectral slope by multiplying it by the power-law 
E^ with C, = 0.55 (see Fig. 7). For calculation of the optical 
thickness ratio we used Ecyc,o = 66.5 keV and velocities on 
the bottom of line-forming region /3o = 0.1 and j3o = 0.4 for 
revolution numbers 1134 and 1138, respectively. During rev¬ 
olution 1134, the luminosity of GX 304-1 is very close to its 
critical value L* « (2 3) x 10^^ erg s“^ (Mushtukov et al. 

2015b), while during revolution number 1138 it is far be¬ 
low the critical luminosity. This fact explains our choice of 
Pq for the considered data sets (see also eq. (8)). The the¬ 
oretical opacities were averaged over the pulse period with 
parameters i = 70° and a = 10° (see eq. (12)). We see that 
theoretical ratio of the opacities shows qualitatively similar 
wave-like structure seen in the data. 

The cyclotron energy i?cyc,o = 66.5 keV is higher than 
line centroid energies measured from the data at any sub- 
critical luminosity states (see Fig. 6 and Tablet). This is 
well explained by our model: we always see a red-shifted 
line absorption feature and only at critical luminosity the 
red-shift (due to the Doppler effect) turns to zero, but the 
estimated line centroid energy might be still red-shifted due 
to the complicated shape of the line. At sufficiently low mass 


accretion rate corresponding to luminosities below 
even the resonant optical thickness drops below unity. At 
this point the cyclotron line caused by scattering in the ac¬ 
cretion flow disappears and the cyclotron line from the NS 
surface at i5 ~ Scyc.o can be observed in the spectrum. We 
thus expect that the positive correlation between the line 
centroid energy and the accretion luminosity discussed in the 
present paper should be replaced by a negative correlation 
at sufficiently low luminosities. Here we do not discuss the 
line formation in a region where the plasma flow is stopped 
by Coulomb collisions and assume that the initial spectrum 
does not contain scattering features. In that sense, the the¬ 
ory does not give a complete model for the line formation, 
however, it points at the effects that have to be taken into 
account and gives a qualitatively satisfactory description of 
the current observational data. 


4 SUMMARY 

We have shown that the cyclotron line absorption feature in 
the spectra of sub-critical XRPs can be produced by scat¬ 
tering off moving electrons in the accretion channel. The 
variations of the line position with the accretion luminos¬ 
ity is naturally explained by changes of the plasma velocity 
profile near the NS surface, which is affected by the radia¬ 
tion flux from the hotspot. The presented model has several 
characteristic features: 

(i) the line centroid energy is positively correlated with 
the luminosity, 

(ii) the line centroid energy is always redshifted from the 
true cyclotron energy corresponding to the surface magnetic 
field at the NS pole, 

(hi) the cyclotron line width is also positively correlated 
with the luminosity, 

(iv) the line is expected to be asymmetric with excess at 
its red wing and it can hardly be approximated by widely 
used Lorentzian or Gaussian profiles, 

(v) the phase-resolved cyclotron line centroid energy and 
the width are expected to be negatively and positively cor¬ 
related with the pulse intensity, respectively. 

In order to get the exact shape of the cyclotron line a self- 
consistent model of the line-forming region including radi¬ 
ation transfer and hydrodynamics is required, which is be¬ 
yond the scope of the present work. Our theory is shown 
to be in good qualitative agreement with the available ob¬ 
servational data on the sub-critical XRP GX 304-1 and has 
profound implications for the interpretation of the data on 
the cyclotron lines observed in XRPs. 
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